Momentum relaxation in a semiconductor proximity-coupled to a disordered s-wave 
superconductor: effect of scattering on topological superconductivity 



Roman M. Lutchyn, 1 Tudor D. Stanescu, 2 and S. Das Sarma 3 

1 Station Q, Microsoft Research, Santa Barbara, CA 93106-6105 
2 Department of Physics, West Virginia University, Morgantown, WV 26506 
3 Condensed Matter Theory Center, Department of Physics, University of Maryland, College Park, MD 20742 

(Dated: compiled October 27, 2011) 

We study the superconducting proximity effect between a conventional semiconductor and a disordered s- 
wave superconductor We calculate the effective momentum relaxation rate in the semiconductor due to pro- 
cesses involving electron tunneling into a disordered superconductor and scattering off impurities. The mag- 
nitude of the effective disorder scattering rate is important for understanding the stability of the topological 
(chiral p-wave) superconducting state that emerges in the semiconductor, since disorder scattering has a detri- 
mental effect and can drive the system into a non-topological state. We find that the effective impurity scattering 
rate involves higher-order tunneling processes and is suppressed due to the destructive quantum interference of 
quasi-particle and quasi-hole trajectories. We show that, despite the fact that both the proximity-induced gap and 
the effective impurity scattering rate depend on interface transparency, there is a large parameter regime where 
the topological superconducting phase is robust against disorder in the superconductor. Thus, we establish that 
the static disorder in the superconductor does not suppress the proximity induced topological superconductivity 
in the semiconductor. 

PACS numbers: 03.67.Lx, 71.10.Pm, 74.45,+c 



Introduction. The possibility of engineering Hamiltonians 
that exploit the properties of the interface between two differ- 
ent materials has recently attracted a lot of attention. There 
are many proposals that exploit magnetic, superconducting 
and other types of properties for spintronics and quantum in- 
formation purposes^. In particular, the prospect of realiz- 
ing exotic topological chiral p-wave superconducting states 
carrying Majorana fermions at the interface between a semi- 
conductor with Rashba spin-orbit coupling and a conventional 
s-wave superconductor in sandwich structures^ 2 -' is very in- 
triguing. The basic concept underlying these proposals is that 
electrons tunneling between different materials inherit their 
physical properties. For example, electrons virtually propa- 
gating in the superconductor "feel" superconducting correla- 
tions. This is the basic idea underlying the superconducting 
proximity effect which is used for realizing the topological 
p + ip superconducting state at the interface. In most of the 
previous studieP^, the s-wave superconductor has been con- 
sidered in the clean limit, and the effect of the superconductor 
disorder on the induced state was not addressed. However, 
it is well-known that impurity scattering in the active sys- 
tem (i.e. in the semiconductor) is detrimental for topological 
superconductivity^^. While semiconductors can be grown 
very clean, most ordinary s-wave superconductors (e.g. Al or 
Nb) are disordered and have very short mean free path I. This 
motivates us to revisit the basics of the superconducting prox- 
imity effect and take into account the effects of disorder in the 
superconductor. As pointed out in Ref . [18| superconducting 
disorder might act similar to impurities in the semiconduc- 
tor, see Fig[lJ>. Thus, the effect of superconducting disorder 
on the stability of the topological phase is an important open 
question which we investigate in this paper. The issue is of 
both conceptual and practical importance: On the conceptual 
side it may appear that disorder residing in the superconductor 
could be detrimental to the semiconductor superconductivity 



since Anderson theorem ruling out the immunity of s-wave su- 
perconductivity to non-magnetic disorder may not necessarily 
extend to topological superconductivity where time-reversal 
invariance may be explicitly broken (e.g. by an external mag- 
netic field creating the Zeeman spin splitting), and on the prac- 
tical side, disorder in the superconductor, if it turns out to be 
detrimental to the proximity-induced superconductivity in the 
semiconductor, may simply completely destroy the topologi- 
cal phase. We note that what is important here is the relative 
magnitude of the proximity-induced superconductivity in the 
semiconductor compared with the proximity-induced momen- 
tum relaxation rate. 

The robustness of topological superconducting phases 
against disorder has been investigated within the simple model 
of one-dimensional spinless p-wave superconductors hosting 
Majorana zero-energy modes at the ends 031 UU and 033 The 
presence or absence of these exponentially localized Majo- 
rana modes defines a topological or non-topological phase. 
The phase boundary between topological and non-topological 
phases is approximately given by Ai n dT ~ 1 where Ai n d and 
r are the proximity-induced gap in this spinless p-wave model 
and the impurity scattering time, respectively. In this paper we 
assume that the semiconductor is clean (we refer the reader 
to Refs. [T3lfT6l IT8l and [19] for more details on how disorder 
in the semiconductor affects topological superconductivity in 
the sandwich structures) and, thus, r is entirely determined 
by the disorder scattering in the superconductor. One can ask 
the following question: is it poss ible to re alize a topological 
phase in the sandwich structured ^ 8 ! 1 U 12 ! 15 ] given that both the 
proximity-induced gap and the disorder scattering rate in the 
semiconductor induced by the superconductor impurities de- 
pend on the interface transparency ? This question is partic- 
ularly significant in view of the topological superconducting 
phase in the semiconductor being equivalent to a spinless p- 
wave superconductor. 



In this paper, we consider a simple model for the semicon- 
ductor/superconductor heterostructure and calculate the mo- 
mentum relaxation rate due to electron tunneling into a dis- 
ordered superconductor. We find that the impurity scatter- 
ing rate r _1 in the superconductor is quite small (i.e. much 
smaller than proximity-induced gap Ai n d). The reason for 
that is two-fold: first, the scattering rate r _1 vanishes in the 
lowest order perturbation theory in tunneling t and involves 
only higher-order processes; second, it is further suppressed 
due to the destructive quantum interference between quasi- 
particle and quasi-hole trajectories by a factor l/pp£, with 
Pf and £ being the Fermi momentum and coherence length 
in the disordered superconductor, respectively. The impor- 
tance of quantum interference effects for higher-order tunnel- 
ing processes has been previously discussed in the literature in 
the context of two-electron tunneling 20 22 . Therefore, we con- 
clude that the condition for the existence of topological super- 
conductivity can be satisfied even in the presence of substan- 
tial superconducting disorder 23 . In the rest of the paper, we 
present a detailed calculation supporting this conclusion. We 
emphasize that, in addition to being important for the realiza- 
tion of topological superconducting phases hosting Majorana 
fermions, our result is also very general and applies to all other 
superconducting heterostructures involving proximity effect. 

Theoretical model. We consider a two-dimensional semi- 
conductor in the proximity to an s-wave superconductor as 
shown in Fig.[T] The Hamiltonian for the semiconductor reads 
(R=l) 
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where S is 2D area of the semiconductor, m* is its effective 
mass, fx is the chemical potential, cxr is the Rashba spin-orbit 
coupling strength, \£(r) = and are Pauli matrices 

acting on the spin degree of freedom. The Zeeman splitting V z 
can be proximity-induced due to the presence of a ferromag- 
netic insulator (not shown in Fig.[l])) in 2D proposals 5 or due 
to an applied in-plane magnetic field V z = QsmI^bB with #sm 
being the ^-factor in the semiconductor, see Refs. EHU and 
[241 for more details on the relevance of Eq. ([T} for generat- 
ing topological superconductivity in generic semiconductor- 
superconductor heterostructures. 

The superconductor can be described by the BCS model 
with Hs being the corresponding mean field Hamiltonian. To 
include disorder effects it is convenient to use exact eigen- 
states formalisrrP3SD. \ n the normal state, single-particle en- 
ergies e n and wavefunctions n (x) in the superconductor are 
defined by the following one-body Schrodinger equation: 
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where V(x) represents a particular realization of the disorder 
potential. We assume here that the terms breaking time re- 
versal symmetry in the superconductor are small, but not in 
the semiconductor, either due to the large difference in the 
^-factors in semiconductor and superconductor 6 8 or because 
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FIG. 1 . (Color online) a) Energy spectrum of the semiconductor with 
Rashba spin-orbit coupling and Zeeman splitting. The position of the 
chemical potential corresponding to the helical state, b) Tunneling 
of electrons into the superconductor leads to a momentum relaxation 
rate in the semiconductor F. 



Zeeman splitting in the semiconductor is proximity-induced 
by ferromagnetic insulator and the amplitude for tunnelin g of 
electrons from superconductor to ferromagnetic insulatoP^ 
is small and can be neglected. In this case, mean field BCS 
Hamiltonian can be diagonalized using the following Bogoli- 
ubov transformation: 



7 t 
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= J d 3 x [U n (^ a (x) - ^ n (x)i a (x)] (3) 



Here the transformation coefficients U n (x) and V n (x) are 
given by the corresponding solution of Bogoliubov-de Gennes 
equation. For spatially homogenous superconducting gap Ao, 
the functions U n (x) and V n (x) can be written as J7 n (x) = 
u n (j) n (x) and V n (x) = v n (j) n (x.). The coherence factors u n 
and v n are given by 



2 V E n 



and vt 
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Here E n = \J e\ + Aq; e n and n (x) are exact eigenval- 
ues and eigenfunctions of the single-particle Hamiltonian ([2]) 
which can be chosen to be real. 

The tunneling Hamiltonian between semiconductor and su- 
perconductor reads 



H t = J2j d3xdV (r(x,r , )^t(x)^(r , )+H.c.) , 



(4) 



where x and r' denote the coordinates in the superconductor 
and semiconductor, respectively, and T(x, r 7 ), in the limit of 
a barrier with low transparency, is defined as 



T(x, y')=U 2 {y-y')8{z) 
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see Refs. l27l and 28 for details. 

Momentum relaxation rate. We now calculate the scatter- 
ing rate of an electron in the initial state \p, a) where p, a are 
electron momentum and spin, respectively, into a state \p' ', a). 
Due to the proximity to the disordered superconductor, mo- 
mentum in the semiconductor is not a good quantum number 
anymore and, as a result, levels \p, a) will have some broaden- 
ing T. Since the superconductor is a good metal, its disorder 
can be well-approximated by short-range impurity scattering. 
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Without any loss of generality, we first study the case when 
ctR and V z are zero and then generalize our results to the heli- 
cal regime at the end of the paper. 

The scattering rate of an electron with momentum p can be 
calculated using Fermi's Golden rule: 



(6) 
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with A p ^ p > being the amplitude for scattering to the state p', 
which can be calculated perturbatively in tunneling t 
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It is easy to show that there is no contribution to the amplitude 
in the lowest order of perturbation theory (see also Ref. 03]) 
and one has to consider higher order processes. The lowest 
non-zero contribution appears in the second order in H t . Af- 
ter simple algebra, one finds that the amplitude in the second 
order in \t\ reads 

A p,p' = 1 J2 y i ^i^W^2^T(x 1 ,ri)T(x 2 ,r / 2 ) (8) 
x(p J ^l^/(n)|0)(0|^//(r 2 )|p / ,(7) 2 n 2 ^(#i)0 TO (x 2 ). 

The factor e n j [e 2 n + Aq) appears due to the cancelation be- 
tween particle and hole contributions to the amplitude A PiP f. 
The amplitude ^ depends on eigenenergies e n and eigen- 
functions (j> n {x) obtained for a particular realization of the 
disorder potential. Therefore, one needs to average the rate 
T over different disorder realizations. The exact eigenstates 
formalism is very convenient here as we will show below. Al- 
ternatively and equivalently, one could do a diagrammatic cal- 
culation of the imaginary part of the self-energy by doing a 
perturbative expansion in tunneling and summing ladder di- 
agrams due to the disorder as shown in Fig. [2] For s-wave 
superconductors these two approaches are equivalent. 
We proceed by first introducing the identity 

(j)n(x l )(l) n (x2) = / d^p^K^, X U X 2 ), 

where K(£,x 1: x 2 ) = Z) n £(f £n)^n(^i)</>n(>2) and 
reducing disorder averaging of the rate F to finding cor- 
relation functions of K(£, a?i, x 2 ). After straightforward 
manipulations, the disorder-averaged scattering rate becomes 
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FIG. 2. (Color online) Diagrammatic perturbation theory in the 
tunneling between semiconductor and superconductor. Disorder- 
averaging is performed at each order in tunneling t. The thick solid 
line represents disorder- averaged Green's function in the supercon- 
ductor G(p, uj). The bottom diagram corresponds to irreducible con- 
tributions as far as disorder averaging is concerned and is calculated 
in the manuscript using exact eigenstates formalism. 



The reducible part can be easily calculated by relating 
(iff (xi, X2)) to the ensemble-averaged normal-state Green's 
function: (iQ(xi,x 2 )) = -£lm(G^(xi,x 2 )) = v F f 12 . 
(Upon averaging over disorder, one can neglect the energy 
dependence of the density of states here, i.e. (^f(O) = 
v F . The function f 12 is given by fi 2 = (e* k ^ Xl_X2 ^)FS 
with (...)fs being the average over electron momentum on 
the Fermi surface. For 3D systems the function f 12 is 
equal to f 12 = sin fe ( ^J x ^ 2|) .) The irreducible part 
(iQ 1 (xi,x 2 )iQ 2 (x3,X4))i r can be expressed in terms of 
the classical diffusion propagators - diffusons and Cooper- 
ons, see, for example, Ref. l29l . Assuming that time- 
reversal symmetry is preserved in the superconductor, diffu- 
sons and Cooperons coincide, 7\,(xi,X2) = V^(xi : x 2 ) — 
7^(xi, x 2 ), and the irreducible part of the correlation func- 
tion ( fTO] ) reads 



(K^ (xi , x 2 )iQ 2 (x 3 , x 4 )) ir 



(11) 



= -^ Re [/l4/23^|^ 1 |(xi,X 3 )+/l 3 /24^|^ 1 |(xi,X 4 )] . 

Using these results, we can now perform integrals over £ and 
£' in Eq. ([9]). Given that the reducible part of the correlation 
function is independent of energy, the energy integrals vanish. 
Thus, irreducible terms do not contribute to the momentum 
relaxation rate, see also diagrammatic calculation in Ref. [TTl 
The contribution of the irreducible part is proportional to 



/ de J d ^ e + Ag e + Ag {K({xi ' X2)K ? (a;3 ' Xa)) - f(x,x,) = J d *'J 



Here the brackets (...) denote averaging over different real- 
izations of the random potential V{x) in the superconductor. 
The correlation function (K^ 1 (x 1 ,X2)K^ 2 (x s ,X4)) consists 
of reducible and irreducible parts, 

(if ?1 (x 1 ,x 2 )^ 2 (x 3) x 4 ))= (10) 
(K ix (xi , x 2 ))(K i2 (x 3 , x 4 )) + (K^ (xi , x 2 )^ 2 (x 3 , x 4 )) ir . 
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where P(t,x, x') is a solution of a 3D diffusion equation: 
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^ /(47rL>t) 3 / 2 . Here D is the dif- 



P(t,x, x') = exp ^ 

fusion constant = lvp/3 with Z and being the mean- 
free path and Fermi velocity, respectively; £ = is the 
effective coherence length in disordered superconductors and 
£o = vp/Ao. Thus, the expression for the rate (r) ^ now 
becomes 



p' ^ i=l.A 

x T(xi, ri)T(x 2 , r 2 )T(x 3 , r^)T(x 4 , r^) 
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Taking into account Eq. ([5]), one can now compute the spatial 
integrals. The integrand is a quickly decaying function and 
converges at short length scales p^ 1 where pp is the Fermi 
momentum in the superconductor. Thus, the dominant contri- 
bution comes from x\ ~ x±, X2 ~ x% and x\ ~ x%, X2 ~ x±, 
respectively, and one finds that (T) = (Ti) + (I^): 

(Ti)=E^-^ ^^/ ^^ WXri ^(ri,r 2 ) 
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Here pf, M is the Fermi momentum in the semiconductor, and 
we have assumed that Pf£, ^ 1. It is convenient to re-write 
the above expressions for the scattering rate (r) in terms of 
the level broadening 7 = 27t^f|^| 2 in the semiconductor due 
to the presence of a bulk metal. Then, one can estimate the 
upper bound on (T) to be 



fa 3 7 2 m* 1 
A m p F £ 



(16) 



As follows from Eq. ([16]), the momentum relaxation rate due 
to the presence of a disordered superconductor is proportional 
to 7 2 whereas the proximity-induced gap is of the order of 7. 
Furthermore, (T) is additionally suppressed due to quantum 
interference effects by a nontrivial (and non-obvious) factor 

i/pf£<i. 

Our results for the momentum relaxation rate can be quali- 
tatively explained as follows. From Eq. fl2] ), one can see that 



the rate (T) is proportional to the probability of a quasiparti- 
cle to return to the junction within the time Aq 1 which is the 
time an unpaired quasiparticle can spend in the superconduc- 
tor in a virtual state. This introduces a length scale in the prob- 
lem £ oc y/D/Ao above which return probability is exponen- 
tially suppressed. Therefore, one can think that the effective 
size of the system relevant for this process is of the order £. 
The momentum relaxation rate generated by a superconduct- 
ing layer of size £ is proportional to the attempt frequency 
Tin ~ \t\ 2 VF times the probability to return to the semicon- 
ductor within time Aq 1 . The latter is given by 7 out /Ao <C 1 
where 7 out ~ \t\ 2 v F D (Here \t\ 2 /^ 3 has dimension E 2 and 
v F D t; 2 scales as l/E). Combining all the terms, we finally 
recover the expression for the rate (T) given by Eq. ( [T6| ). In 
particular, these arguments explain why the rate is suppressed 
by a factor 1/pfC 

We now take into account the effect of Rashba spin-orbit 
coupling and Zeeman splitting in the semiconductor and dis- 
cuss momentum relaxation rate in the helical phase. The situ- 
ation we are interested in is when the chemical potential is 
in the gap as shown in Fig. [ij. Thus, the momentum re- 
laxation rate is determined by the scattering amplitude be- 
tween different momenta on the helical Fermi surface which 
is given by A PiP > = - (p, -| H t ^H t \p', -). Here \p, ±) 
denotes the state of an electron on the Fermi surface with 
a particular chirality ±. The results of the calculation can 
be straightforwardly obtained by repeating the steps above 
but are not particular illuminating. For experimentally rele- 
vant parameters, the rate (T) has to be multiplied by a fac- 
tor 0(1) which originates from the change in the matrix el- 
ements, as well as the density of states. Therefore, for prac- 
tical purposes one can use Eq. [16] in 2D system. It is worth 
pointing out, however, that in the limit m*a 2 ^> V z there is 
an additional suppression of the impurity scattering rate be- 
cause the Berry phase of the Fermi surface is equal to 7r up 
to corrections V z /m*a 2 <C 1. This suppression of the elastic 
backscattering is particularly important for one-dimensional 
helical nanowires. Overall, the rate fT6] ) should be considered 
as an upper bound on the effective impurity scattering rate due 
to superconducting disorder. We can now estimate (r)/A in d. 
Taking InAs as a semiconductor and Al (or Nb) as a super- 
conductor, we find that 1/pf£ ~ 10 -3 and m*/m « 0.04. 
The effective proximity-induced gap varies from Ai nc j ~ 7 
for m*a 2 > V z and A ~ j^m*a 2 /V z for m*a 2 <C V z . 
Assuming that 7 ~ A and taking the pessimistic numbers 
for spin-orbit coupling m*a 2 /V z ~ 0.1, one finds that the 
ratio of the scattering rate (T) fT6] ) to the proximity-induced 
gap Ai n d ~ 7 is small (r)/Ai n d ~ 10 -2 . This ratio involving 
the upper bound of the momentum relaxation rate is thus very 
small and we, therefore, conclude that the topological super- 
conducting phase emerging at the interface is robust against 
disorder in the superconductor. 
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